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SECTION  1 . INTRODUCTION 

Strong  shocks  in  air  reflect  off  rigid  surfaces  at  many 
times  (5-13  times)  their  original  pressure.  Opposing  shocks 
of  equal  strength  behave  just  as  a reflected  shock , leading 
to  very  high  pressures  at  the  initial  point  of  contact.  If 
two  shocks  of  100  psi  collide,  the  resulting  peak  pressure  is 
around  500.  For  two  1000-psi  shocks,  the  peak  pressure  jumps 
to  8500  psi.  This  suggests  that  the  area  covered  by  a given 
overpressure  from  two  simultaneous  blast  waves  may  be  consid- 
erably large::  than  twice  the  area  of  a single  burst. 

Yet,  while  strong  shock  reflection  factors  are  impressive, 
there  is  reason  to  believe  that  such  high  values  do  not  extend 
far  beyond  the  initial  point  of  contact,  and,  further,  two 
unequal  shocks  may  interact  in  a much  less  impressive  way. 

Strong  blast  waves  are  extremely  transitory,  and  pressures, 
densities,  and  flow  rates  behind  each  blast  front  drop  off 
exceedingly  rapidly. 

Careful  considerations  of  such  blast  wave  interactions 
lead  directly  to  three-dimensional  geometries  which  very  much 
inhibit  the  accuracy  and  practical  resolution  achievable  with 
canonical  numerical  methods.  The  following  series  of  estimates, 
without  benefit  of  rigorous  modeling  and  detailed  numerical 
calculations,  are  intended  to  bound  the  expectations  for  enhanced 
coverage  by  means  of  simultaneous  blasts.  Cotiparison  is  also 
made  with  the  LAMB  procedure,  as  applied  to  two  simultaneous 
bursts  and  the  overpressures  along  the  line  joining  their  centers. 

Low-Altitude  Multiple  Burst. 


SECTION  2.  SEVERAL  ESTIMATES 


The  peak  overpressure  from  a single  burst  on  the  surface 
is  well  approximated  with  the  formula  '1] 


3300  W . 192  rfi 

4p  - — 3-  + -pTj-  psi 


(l) 


with  W the  yield  in  megatons  and  R the  distance  in  kilofeet. 


A normally  reflected  shock  reaches  a reflected  pressure 
enhanced  by  the  stagnation  of  the  flow,  and  results  in  pres- 
sures for  a strong  shock  much  more  than  double  the  incident 
shock  pressure.  An  approximate  shock  reflection  factor  for 
an  ideal  gas  of  specific  heat  ratio  y is  given  by 

Apr  4yPQ  + (3y  - 1)AP 

AP“~  R = 2yP  + (y  - 1) AP  (2) 

o 

where  AP  is  the  incident  overpressure,  Apr  is  the  reflected 
overpressure,  and  PQ  is  the  ambient  pressure  (14.7  psi) . For 
sea  level  air,  1.1  < y < 1.7  (y  = 1.3  for  AP  = 1000  psi) (1]. 


A more  exact  fit  (within  3 percent)  to  this  reflection 
factor  is  provided  by  the  formula  below,  which  accounts  for 
the  nonideal  gas  properties  of  sea  level  air  [2]. 


W — — 

This  formula  agrees  to  within  10  percent  with  the  accepted 
average  of  the  nuclear  test  data  which  in  turn  are  90  percent 
contained  by  a spread  of  ±50  percent  at  pressures  above  40  psi. 


.TW-'jyt* 


RF  = 


0.002655AP 


-9  2 + 2 

1 + 0.0001728AP  + 1.921  x 10  *AP* 

0.004218  4-  0 . 04834AP  + 6.856  x 10"6AP2 
1 + 0. 007997AP  + 3.844  x 10"6AP2 
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At  the  point  where  two  simultaneous  spherical  blast  waves 
meet,  the  peak  overpressure  can  be  fairly  precisely  described 
with  these  two  approximations  (Equations  1 and  3) . For 
example  (as  in  Figure  1) , for  two  1-MT  simultaneous  surface 
explosions  separated  by  6860  ft,  the  blast  waves  meet  when 
each  has  a peak  overpressure  of  116  psi  (Equation  1) , and 
these  shocks  reflect  at  the  point  of  contact  to  a pressure  of 
600  psi  (a  reflection  factor  of  5.18) (Equation  3).  The  value 
of  600  psi  from  a single  1-MT  surface  burst  occurs  at  1840  ft. 
For  nonsimultaneous  bursts,  a separaticn  of  3680  ft  would 
cover  a line  target  with  more  than  600  psi  everywhere.  The 
separation  distance  at  which  the  interacting  shocks  reflect 
to  600  psi  (6S6C  ft)  is  nearly  twice  as  long. 

The  area  associated  with  the  region  of  enhanced  blast 
pressures  is  a thin  lens  about  the  point  of  contact  between 
spherical  (or  hemispherical)  shocks,  and  is  only  a small 

* 

fraction  of  the  area  covered  by  the  individual  blast  wr-es. 

However,  the  use  of  the  larger  distance  (6860  ft) , where 
the  shocks  reflect  to  just  600  psi  as  an  effective  kill  distance 
for  a line  target,  is  quite  incorrect.  The  geometry  of  the 
situation  (Figure  2)  suggests  that  as  the  two  shocks  pass 

It  is  very  wrong  to  assume  that  if  the  separation  distance 
were  increased  by  a factor  of  two,  the  area  coverage  would  be 
correspondingly  increased  by  a factor  of  four. 


, througn  ^uch  other-  they  continue  to  expand  spherically,  and 

* 

; so  drop  rapid. _ in  pressure  from  that  peak  value  of  first  con- 

• tact  l).  if  their  peak  interaction  gives  only  600  psi, 

then  all  iiw  region  between  burst  points  that  lies  beyond  a 

; radius  from  each  burst  point  of  1840  ft  would  see  less  than 
600  (Figure  3).  That  is,  as  in  Figure  3,  a total  of 
3200  ft  would  experience  two  shocks — both  less  than  600  psi. 

Of  create/"  interest  is  the  separation  distance  that  exposes  jj 

the  fer.t  '• .*  * :ine  target  to  more  than  600  psi  (or  some  other  peak  | 

overpressure  of  interest)  . For  this  we  need  to  know  how  rapidly  i 

| these  diverging  shocks  decrease  in  overpressure  as  they  expand  I 

• into  each  other.  Since,  on  first  contact  with  the  opposing  ji 

V ) 

’ blast  w-/ves  the  transmitted  shock  starts  at  the  peak  reflection  j 

l 

\ v/'.lue,  the  initial  value  It  well  defined  (Equation  3),  but  there  j 

j r’-'es  not  exist  an  equally  simple  or  direct  formula  for  predicting 

i the  subsequent  decay  rate  as  the  two  blasts  penetrate  each  other. 

! | 

• Numerical  blast  calculations  provide  detailed  descriptions 

! ; , 

i of  the  pressure  field  into  which  each  transmitted  shock  is  i | 

i 1 

i expanding  and  how  it  behaves  in  both  space  and  time  [3,4,5].  j 

| For  most  applications  to  hard  targets  (e.g.,  for  600-psi  ' 

' trenches) , a simple  strong  shock  model  would  suffice.  In  any  ; 

case,  while  peak  pressure  decays  with  distance  (as  in  Equation  1) 
approximately  in  proportion  to  the  inverse  cube  of  the  distance, 
the  interior  of  the  fireball/blast-wave  follows  with  a pressure 
about  one- third  to  one-half  of  the  shock  value  (Figure  4) . In  | 

fact,  both  peak  pressure  and  interior  pressure  decrease  in  pro- 
portion to  the  inverse  of  the  time  measured  from  the  instant 

★ 

, of  explosion.  An  approximate  relation  is  given  by  [2] 


1 


Two  simpler  but  less  exact  forms  for  the  pressure-time  rela- 
tion  in  a nuclear  blast  wave  may  be  found  in  Reference  1, 

pp.  180-181. 
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AP ( t,T) 


( 148000  \ 

\Q. 135  + t) 


0.417  + 0.583 


(4) 


in  which  t is  the  time  in  msec  (after  burst) , T is  the  shock 
arrival  time  (msec)  for  1 MT,  (t  - T) , D is  the  positive  phase 
duration  (msec)  (the  time  during  which  the  blast  pressure  is 
greater  than  ambient),  and  f(t)  is  an  empirical  adjustment  fit 
of  secondary  importance. 


_ 100  + 6. 72t  + 0.00581t2 
fit;  - y— 

100  + 18. 8t  + 0.0216tz 


(5) 


The  essential  time  behavior  of  the  pressure  is  illustrated 

by  this  approximation:  The  dominant  behavior  is  a decay  almost 

-1  15 

linear  in  time  (more  precisely  as  ^ t * ) with  a very  sharp 
drop  just  behind  the  shock  front  to  a value  of  about  40  percent 
of  that  at  the  front  {'v  t ) . 

One  possible  approximation  is  to  assume  that  the  transmitted 
shock  continues  to  generate  the  same  peak  reflection  factor  as 
it  continues  to  expand  and  decrease  inside  the  other  blast  wave. 
This  is  likely  a gross  overestimate  of  the  off-peak  pressures, 
since  one  might  better  use  reflection  factors  appropriate  to 
the  transmitted  shock  as  it  continues  to  expand  and  decrease. 

The  original  blast  that  it  is  running  into  also  continues  to 
expand  and  decrease  in  pressure.  A more  correct  approximation 
will  account  for  this  double  decay,  but,  for  the  moment,  consider 
several  simple  approximations  for  the  transmitted  peak  over- 
pressure . 

Figure  5 shows  several  choicee  for  these  transmitted 
pressures  in  the  particular  case  of  two  simultaneous  surface 
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bursts  (1-MT)  separated  by  4500  ft.  Table  2 illustrates  the 
development  of  numerical  values  for  each  case. 

(1)  The  peak  interaction  is  defined  as  the  local 

pressure  in  front  of  the  transmitted  blast  wave 
multiplied  by  the  initial  reflection  factor  for 


the  colliding  shock  fronts  (reflection  factor 
RF  - 7.0) . 

(2)  The  peak  interaction  pressure  is  defined  as  the 
local  pressure  in  front  of  the  transmitted  blast 
wave  multiplied  by  the  reflection  factor  appropriate 
for  a shock  of  strength  equal  to  that  of  the  second 
shock  if  it  had  expanded  (unreflected)  to  that  dis~ 
tance  beyond  the  initial  contact  point. 

(3)  The  peak  interaction  pressure  is  defined  as  the 
local  pressure  in  front  of  the  transmitted  blast 
wave  multiplied  by  the  reflection  factor  for  a 
shock  of  peak  pressure  equal  to  that  pressure. 

(4)  The  peak  interaction  pressure  is  defined  as  the 
unreflected  (incident)  shock  pressure  multiplied 
by  the  reflection  factor  appropriate  to  the  time 
of  first  contact  (RF  * 7.0). 


(5) 


The  peak  interaction  pressure  is  defined  as  the 
unreflected  (incident)  shock  pressure  multiplied 
by  the  reflection  factor  appropriate  for  that 
reduced  shock  strength  as  it  expands  in  an  undis- 
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Table  2.  1-MT  Simultaneous  Burst  Interaction  (Separated  by  4500  ft)  Approximations 

to  Peak  Overpressure 


mwwi 


None  of  these  approximations  account  for  the  expansion 
of  both  shocks  as  they  pass  through  each  other.  The  most 
reasonable  approximations  are  (3)  and  (2) , but  the  others 
are  shown  for  comparison  (Figure  5) . 

The  extra  line  target  coverage  for  this  example  (1-MT 
bursts  separated  by  4500  ft)  is  listed  in  Table  3. 


Table  3.  Comparison  of  Range  and  Area  Coverage  Increases 
for  Five  Different  Approximations  to  the  Interaction 
of  Two  Simultaneous  1-MT  Surface  Bursts  Separated  by  4500  ft 


APPROXIMATION 

NUMBER 

PEAK 

OVERPRESSURE 

COVERED 

(psi) 

SINGLE 

BURST 

RANGE 

(ft) 

DISTANCE 
ADDED  BY 
INTERACTION 
(ft) 

3a 

540 

1915 

670 

17 

2a 

610 

1830 

840 

23 

1 

640 

1800 

910 

25 

5 

990 

1540 

1410 

46 

4 

1120 

1470 

1550 

53 

aM0ST  PLAUSIBLE  APPROXIMATIONS. 

The  last  two  are  clearly  overestimates  of  the  effect, 
but  the  first  three  may  also  predict  more  enhancement  than  is 
real. 

Further  estimates  are  given  in  the  following  sections. 
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SECTION  3.  INTERACTION  OF  UNEQUAL  SHOCKS 
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The  preceding  discussion  has  dealt  with  the  case  of  two 
simultaneous  bursts  (equal  strength  shocks)  interacting. 
Interactions  between  shocks  of  unequal  strength  are  also  of 
interest,  since  timing  of  multiple  bursts  may  not  be  exact 
and  the  first  blast  may  be  considerably  weaker  by  the  time 
the  second  blast  meets  it. 


Suppose  two  shocks  meet  when  their  respective  peak  over- 
pressures are  AP^  and  AP2,  and  suppose  the  first  is  stronger 
than  the  second  (AP1  > AP2) . These  shocks,  on  interacting, 
lead  to  an  overpressure  APR,  a density  pR*  and  a resultant 
particle  velocity  UR.  Figure  6 identifies  the  nomenclature 
in  which  two  initial  shocks  have  met  and  have  resulted  in  two 
transmitted  shocks  which  have  velocities  UR1  and  UR2  directed 
away  from  each  other,  and  a particle  velocity  UR  which  is 
positive  to  the  right  if  AP1  > AP2. 


! 


I 

! 

i 

! 

! 


i 
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Figure  6.  Shock  Interaction  Notation  for  Unequal  Shocks 




Actually,  two  states  of  density  and  temperature  separated  by 
a contact  discontinuity  exist  in  the  shock  interaction  region 
(shaded  area,  Figure  6) , but  pressure  and  velocity  are  the 
same  in  both  states,  viz.,  APR,  nR. 
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From  standard  shock  relations/  U^,  and  P2  are 

expressible  in  terms  of  AP^  and  and  are  considered  as 
known  here. 


“i(<1r1>4Pi + *po)  uy^0 

4P2(<l±i>Ap2  ♦ YP0)  in/s^ 

> [(<xrjK + >po]f(<xr)  AP1  + *.]  1 
, [(<xirK + ,po]  [(<xr)  AP2  + vpJ'1  . 


Expressing  the  conservation  of  mass  flow  rate  through  each 
shock,  one  may  write 


PR1<°R1  + UR>  * pl(Ul  + °R1>  ' 


PR2<UR2  - V * P2(U2  + UR2>  * (7> 

The  momentum  change  related  to  the  pressure  jumps  at  the 
shocks  can  be  expressed  as 

APr  - APj  - P10l(0x  + URl)  - PR1UR(0R  + UR1) 

APr  - APj  - P2U2(Uj  ♦ UR2)  + PR20R(UR2  - UR)  . (8) 
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Eliminating  UR1  and  UD,  from  Equations  6 and  7 , one 


can  write 


R2 


Pi  PRi  2 

4PR  ■ AP1  + <U1  - UR> 


and 


^2PR2  2 

fiPR  * AP2  + <°2  + UR» 


(9) 


The  initial  energy  density  (E)  jump  conditions  at  each 
transmitted  shock  require 


and 


ER  - E1  - 


Er'E2 


R+Pl  /l_  _1_\ 

2 ^1  PRl) 

^'^2) 


PR+P2 


(10) 


If  one  assumes  an  ideal  gas  with  ratio  of  specific  heats  (y) , 
then  E *P/[(y-l)pJ.  Solving  for  pR  leads  to 


>R1  _ (y+1)  APR  + (y-DAPj^  + 2yP0 

57"  rv-i(ipR  + 


and 


PR2  _ <Y+1)APr  + <y-l)AP2  + 2yP 
P2  * (y-l)APR  + (y+l)Af;2  + 2yP 


(ID 


Using  Equations  10  and  the  Hugoniot  relations  to  express 
P^t  P')*  and  U2  in  terms  of  AP^  and  AP2  (Equation  6),  one 


can  derive 


APR  « A ± A + C 


APR  * D ± + P 


in  which 


A = Api  + lii  p1<0l  - oR) 


5 -api2  + iPi  + *po)  Pi10!  - V 


iP2  4 ^ P2(U2  + V 


-Ap22  + AP2  + yO  P2(U2  + UR) 


These  equations  can  be  solved  for  APr  by  iterative 
selection  of  the  velocity  UR  (all  other  values  being  pre- 
determined by  the  values  of  YPqp  AP^,  and  AP2) . 

A set  of  example  values  of  the  resultant  peak  overpressure 
APr  for  two  shocks  meeting  when  one  of  them  is  100  psi  is 
illustrated  in  Figure  7 for  two  values  of  the  specific  heat 
ratio.  For  example,  a 100-psi  and  a 500-psi  shock  meet  to 
give  more  than  1500  psi  for  y * 1.4  and  about  2000  psi  for 
y - 1.2. 


Figure  8 illustrates  the  peak  pressure  amplification 
from  the  meeting  of  two  shocks  with  a plot  of  the  ratio  of 
the  resulting  peak  overpressure  to  the  sum  of  the  two  incident 
peak  pressures.  The  amplification  approaches  unity  if  one  of 


21 


the  shocks  is  weak#  and  seems  to  reach  a peak  when  the  ratio 
of  the  two  incident  overpressures  lies  between  3 and  5. 

This  solution  for  the  resultant  peak  pressure  from  the 
collision  of  two  unequal  shocks  suggests  yet  another  approxi- 
mation to  the  peak  overpressure  distribution  for  our  example 
of  two  simultaneous  1-MT  bursts  separated  by  4500  ft  (Figure  5) . 
The  interior  overpressure  (behind  the  shock  front)  that  a 
transmitted  shock  encounters  can  be  assumed  to  act  similarly 
to  a shock  of  that  overpressure,  and  the  resultant  overpressure 
computed  (by  means  of  Equations  12  and  13) . This  approximation 
is  tabulated  in  Table  4 and  illustrated  in  Figure  9 (labeled 
Case  6) . 

This  approximation,  using  unequal  shock  properties,  is 
not  very  different  from  the  simple  average  of  the  peak  over- 
pressures from  the  earlier  bounding  cases— Cases  3 and  5— 
which  use  the  interior  blast  wave  pressure  multiplied  by  the 


normal  shock  reflection  factor  for  that  pressure  (Case  3)  and 
the  single  shock  peak  overpressure  multiplied  by  its  normal 
reflection  factor  (Case  5) . This  average  is  also  shown  in 
Figure  9 as  Case  7. 
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SECTION  4.  COMPARISON  WITH  THE  LAMB  PROCEDURE 

A more  careful  accounting  for  the  usual  conservation  of 
mass,  momentum,  and  energy  during  these  shock  interactions 
requires  some  further  assumptions  and  more  geometry  and  arith- 
metic. Such  an  attempt  is  the  basis  of  extensive  calculations 
and  predictions  at  the  Air  Force  Weapons  Laboratory  with  a 
program  referred  to  as  LAMB  [6]. 

The  essential  features  of  this  procedure  are  as  follows: 
Mass  conservation;  assume: 

- n 

P * P,>  + £ Ap,,  (P  - 0*5  P_)  (14) 

° i=l  1 ° 

in  which  p is  the  air  density,  pQ  the  ambient  (pre-shock)  air 
density,  and  Api  is  the  over-density  (p.^  - pQ)  in  each  blast 
wave . 


Recognize  that  densities  in  strong  shock  fronts  may  rise 
to  more  than  ten  times  the  ambient  density,  but  may  also  fall 
to  less  than  one-tenth  the  ambient  value  inside  the  fireball. 
This  prescription  predicts  the  peak  density  of  two  equal 
colliding  strong  shocks  to  be  only  about  half  of  the  correct 
Hugoniot  value.  (See  Figure  10.) 

The  program,  originally  designed  to  approximate  over- 
lapping bursts  at  altitudes  between  10  and  30  kft  (in  a 
missile  defense  role) quite  arbitrarily  restricts  the  under- 
densities (negative  overdensities)  to  half  (or  more)  of  the 
ambient  density.  Such  a restriction  provides  a reasonable, 
if  unsatisfactorily  empirical  accounting  of  other  disoersive 


irlson  of  Reflection  Factor  (A°R/A°S)  for  Normal  Shocks  and  by  the 
, Method;  Comparison  of  Reflected  Density  and  Density  by  Lambs  Method 


effects  likely  to  occur  in  the  high  temperature  (low  density) 
interior  of  a nuclear  fireball,  but  it  cannot  be  justified 
rigorously  (from  first  principles).  Actually,  some  such  limit 
is  needed  to  prevent  an  undesired  increase  in  net  kinetic 
energy  for  a shock  in  the  very  low  density  interior  of  a strong 
blast  wave  (fireball)  that  would  result  otherwise  from  this 
prescription.  Recognize  also  that  this  expression  is  not  one 
of  mass  conservation,  but  prescribes  density,  or  mass-per-unit 
volume.  It  has  long  been  a favorite  piece  of  magic  for  simpli- 
fied blast  solutions  (many  published  as  serious  contributions) 
to  make  seemingly  innocuous  assumptions  about  density  distribu- 
tions and  then  proceed  to  unfold  a marvelously  consistent 
picture  of  some  blast  wave.  The  rabbit  is  always  already  in 
the  hat  here,  however,  since  density  distributions  are,  in 
fact,  integral  representations  of  the  entire  movement  history 
of  the  blast,  and  the  movements  so  described  are  a direct 
consequence  of  the  acceleration  or  force  (pressure  gradient) 
history.  So  any  density  profile  contains  the  blast  wave 
history  to  that  instant,  and, is  not  a trivially  adjustable 
parameter  of  small  consequence. 

Conservation  of  momentum;  assume: 

. n 

pv  - £ p.v.  (15) 

i-1  1 1 

where  p.  is  the  density  in  the  i blast  wave,  v.  is  the 
particle  velocity  of  the  i blast  wave,  and  p and  v are  the 
resultant  density  and  particle  velocity  in  the  blast  interac- 
tions. This  represents  a local  momentum  density  vector  sum 
without  consideration  of  pressure  impulse  contributions,  which 
is  not  quite  consistent  with  the  usual  assumptions  of  inviscid 
gas  dynamics  for  blast  wave  characterization. 
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Conservation  of  energy;  assume: 


AP  = p . Po  = £ APi  + i ft  P1vi2  - pv2j  (16) 

in  which  AP^  = - PQ  is  the  overpressure  in  the  i blast, 

and  AP  is  the  resulting  overpressure  in  the  interacting  blasts. 

This  prescription  purports  to  convert  excess  kinetic 
energy  from  the  opposing  flows  into  pressure  energy,  a proce- 
dure consistent  with  normal  hydrodynamic  flow  characterizations. 

2 

However,  it  is  inexact  as  an  energy  equivalent  to  add  1/2  pv 
terms  to  overpressure.  The  dimensions  are  correct,  but  the 
compressibility  factor  l/(y-l)  is  missing.  As  a consequence, 
the  reflection  factor  for  two  equal  shocks  resulting  from 
this  prescription  is  slightly  in  error.  For  an  ideal  gas: 

(AP  \ 2yAP  + 4yPQ 

“"/lamb  * Rlamb  ■ (y-I)AP  + 2ypo  • (17) 

More  rigorously  (as  in  Equation  3) , 

AP  (3y-l)AP  + 4yP 

AP~  R “ (y-l)AP  + SyPQ  * (18) 

4 

These  two  formulae  are  compared  in  Figure  10  for  1 < AP  < 10  , 
showing  the  LAMB , procedure  to  be  low  by  a factor  of  7/8  at 
high  overpressures  (y  ■ 1.4). 

The  LAMB  procedure,  when  applied  to  the  previous  example 
(two  simultaneous  1-MT  surface  bursts  4500  ft  apart) , gives 
even  less  coverage  with  high  pressure  than  the  previous 


estimates  using  peak  reflection  factors.  As  mentioned,  the 
LAMB  method  misses  the  peak  pressure  by  a factor  approaching 
7/8,  and  drops  away  from  that  peak  value  very  rapidly  (Figure  11). 
Table  5 compares  the  range  enhancements  from  the  approximations 
previously  discussed  with  that  for  the  LAMB  model.  The  curves 
in  Figure  11  and  the  basic  data  m Tables  2 and  4 are  derived 
from  Reference  1,  but  any  of  the  descriptions  in  References  3, 

4 or  5 would  serve  as  well. 

The  added  line  coverage  for  this  example  with  the  LAMB 
procedure  is  about  11  percent,  a smaller  effect  than  any  of 
the  approximations  given  in  the  previous  section.  Good  agree- 
ment with  experiments  is  claimed  for  the  LAMB  model  when  applied 
to  HE  tests  or  to  the  shock  reflections  on  a nuclear  test  such 
as  PLUMBBOB-PRISCILLA.  However,  the  arbitrary  assumptions  in 
the  model  are  neither  intuitive  nor  physically  correct,  and 
their  effect  on  results  far  from  the  point  of  initial  shock 
contact  remains  unclear. 
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SECTION  5.  SUMMARY  AND  CONCLUSIONS 


S 


| A reasonable  lower  limit  approximation  to  the  combined 

| overpressure  would  appear  to  be  to  multiply  the  pressure  at 

f 

i a point  inside  a single  blast  and  in  front  of  the  transmitted 

shock  from  a second  burst  by  the  reflection  factor  for  a shock 
of  that  interior  pressure  {Case  3) . A similarly  simple  upper 
bound  should  be  provided  by  multiplying  the  single  blast  peak 
overpressure  by  the  normal  reflection  factor  for  that  pressure 
at  distances  beyond  the  point  of  first  contact  for  the  two 
shocks  (Case  5).  A further  procedure , almost  as  simple,  pro- 
viding an  intermediate  value,  uses  the  pressure  predicted  for 
the  interaction  of  unequal  shocks,  based  on  the  pressure  ahead 
of  and  behind  the  transmitted  shock  (Case  6) . 

The  LAMB  procedure  is  not  rigorously  correct,  but  it 
provides  a more  general  and  a more  detailed  treatment  of 
interacting  shocks.  Unfortunately,  it  is  not  clear  whether 
it  overestimates  or  underestimates  the  resulting  pressures, 
although  all  the  predictions  for  this  example  lie  above  the 
LAMB-derived  values.  Accepting  all  the  apparent  uncertainty 
in  these  approximations  to  multiple  shock  interactions,  one 
can  still  conclude  that: 


*•{ 


• The  region  of  enhanced  overpressures  from  two 
simultaneous  separate  blast  waves  is  a small 
fraction  of  the  area  covered  by  each  individual 
blast.  r 


• The  extra  coverage  of  a line  target  with  high 
overpressures  (AP  > 300  psi)  is  of  the  order  of 

o 

10-50  percent,  and, .more  probably,  less  than 
30  percent. 
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The  LAMB  procedure  may#  in  some  cases#  underpredict 
the  peak  pressures  for  interacting  blast  waves. 
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ATTN:  Sec.  Officer  for  J.  Wlsotski 

Pacifica  Technology 

l 

f 

EC&G  Washington  Analytical  Services  Center,  Inc. 

ATTN:  R.  Bjork 

] 

ATTN:  Library 

ATTN:  C.  Kent 

Card,  Inc. 

Physics  International  Co. 

! 

ATTN:  G.  Neidhardt  : 

ATTN:  C.  Vincent 

ATTN:  F.  Sauer 

! 

General  Electric  Co. 

ATTN:  Technical  Library 

Space  Division 

ATTN:  D.  Orphal 

‘ 

ATTN:  M.  Bortrer 

ATTN:  E.  Moore 

ATTN:  L.  Behrmann 

i 

i 

J 

Generel  Electric  Co. -TEMPO 

Center  for  Advenced  Studies 

R 4 D Associates 

>' 

i 

ATTN:  DAS I AC 

ATTN:  R.  Port 

ATTN:  W.  Wright,  Jr. 

* 

J 

} 

J 

Ceneral  Research  Corp. 

ATTN:  C.  MacDonald 

t 

J, 

Santa  Barbara  Division 

ATTN:  J.  Carpenter 

j 

s' 

ATTN:  B.  Alexander 

ATTN:  J.  Lewis 

ATTN:  Technical  Information  Cantor 

* 

i. 

1IT  Research  Institute 

ATTN:  C.  Knowles 

1 

S 

ATTN:  R.  Welch 

ATTN:  A.  Kuhl 

i 

ATTN:  M.  Johnson 

10  cy  ATTN:  A.  Lttter 

1 

ATTN:  Documents  library 

10  cy  ATTN:  H.  Brode 

\ 

? 

Information  Sciences,  Inc. 

R & D Associates 

I 

ATTNl  W.  Dudziak 

ATTN:  H.  Cooper 

1 

Inatitute  for  Defense  Analyses 

Rand  Corp. 

1 

ATTN:  Classified  Library 

ATTN:  0,  Mow 

J.  H.  Wiggins  Co.,  Inc. 

Science  Applications,  Inc. 

ATTN:  J.  Collins 

ATTN:  Technical  Library 

! 

( 

Kaaan  AvlDyne 

Science  Applications,  Inc. 

i 

Division  of  Kaaan  Sciences  Corp. 

ATTN:  J.  Dishon 

ATTN:  E.  Crlaclone 

ATTN:  Library 

Science  Applications,  Inc, 

\ 

| 

ATTN:  N.  Hobbs 

ATTN:  M.  Knssel 

ATTN:  R.  Slavers 

i 

Roman  Sciences  Corp. 

ATTN:  B.  Chambers 

J ^ 

ATTN:  F.  Shelton 

ATTN:  Library 

Science  Applications,  Inc. 

! , 

ATTN:  D.  Maxwell 

i i 

Lockheed  Missiles  & Space  Co.,  Inc, 

ATTN:  D.  Bernstein 

ATTN:  TIC,  Library 

Southwest  Research  Institute 

i 

j 

‘ 

Lovelace  Biomedical  & Environmental  Research 

ATTN:  W.  Baker 

j ; 

Institute,  Inc. 

ATTN:  A.  Wenzel 

| 

ATTN:  R.  Jones 

SRI  International 

j | 

Martin  Marietta  Corp. 

ATTN:  C.  Abrahaoeon 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Terra  Tek,  Inc. 

ATTN:  Library 
ATTN:  A.  Jonea 
ATTN:  S.  Green 

Tetra  Tech.,  Inc. 

ATTN:  Library 
ATTN:  L.  Hwang 

TRW  Oefenae  A Space  Sya.  Group 
ATTN:  R.  Plebuch 
ATTN:  P.  Bhutta 

ATTN:  Technical  Infonuclor.  Center 
ATTN:  D.  Baer 
ATTN:  I.  AXber 
2 cy  ATTN:  P.  Dai 

TRW  Defenae  A Space  Sya.  Group 
San  Bernardino  Operations 
ATTN:  E.  Wong 
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Universal  Analytics,  Inc. 

ATTN:  E.  Field 

Eric  H.  Wang 

Civil  Engineering  Rach.  Fac. 

ATTN:  N.  Bauo 
ATTN:  L.  Bickle 

Weld linger  Assoc.,  Consulting  Englneere 
ATTN:  M.  Baron 

Weldlinger  Assoc.,  Consulting  Engineers 
ATTN:  J.  Isenberg 

Weatlnghouae  Electric  Corp. 

Marine  Division 

ATTN : W.  Volz 
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